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kaon is presented within the framework of light-cone QCD, with perturbative fcr-factorization in- 
cluding Sudakov suppression. The effects of power suppressed sub-leading twists and the gen- 
uine "soft" QCD corrections turn out to be dominant at low- and moderate-energies/momentum 
transfers Q^ . Our predictions agree weh with the available moderate-energy experimental data, 
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I. INTRODUCTION 

Electromagnetic form factors are interesting physical observables in hadronic physics 
which directly provide insights into the hadronic constituents, charge distributions, currents, 
color and flavor within the hadrons. Their precise knowledge is of fundamental importance 
for a realistic and accurate description of exclusive nuclear reactions that serve as ideal 
testing grounds for understanding the dynamics of confinement in QCD that have been 
grappling with physicists ever since the discovery of asymptotic freedom. 

In the last few decades, there has been significant experimental efforts in extracting 
hadronic form factors (e.g., see. (l|-l3| for the charged pion form factors ) from various exclu- 
sive processes. However, in case of the charged kaon form factors their behavior were very 
severely constrained due to absence of reliable experimental data. Since the mid-90s kaon 
photo-/electro-production experiments on reactions such as A{'y, K)YB and A(e, e'K)YB 
(target A, produced hyperon Y and recoil B) have invited some renewed interest in the study 
of kaon form factors, although the existing data is still too limited, restricted only to the very 
low space-like region, as low as —q^ < 0.2 GeV^ [4]. In the time-like region, there are more 
scattered data up to several GeV^ (albeit with very large error-bars) for time-like processes 
such as 7* — > K^K^, extracted from annihilation reactions such J/ip — )■ e"'"e~ — )■ h^h~ 
{h = IT, K, ...) by applying suitable experimental cuts. A compilation of previously extracted 
kaon form factors for q^ = Q^ < 10 GeV^ is given in |2|. Recently, high precision measure- 
ments by the CLEO collaboration with first ever identified time-like kaons for Q^ > 4 GeV^ 
have yielded the following resuhs: |Gi^(13.48 GeV^)] = 0.063 ±0.004(stat)±0.001(syst) and 
Q2 |G'K(13.48GeV2)| = 0.85 ± 0.05(stat) ±0.02(syst) GeV^ [3,]. Note that in this paper, 
we shall use the symbol Gk for the time-like kaon form factor to distinguish it from the 
space-like counterpart Fk- 

Notwithstanding the aforementioned problem of paucity of quality statistics of the exist- 
ing kaon data, the purpose of the paper is an effort to make a prediction for the charged 
kaon form factors using the framework of perturbative factorization p, la]. In this way, 
we hope to throw some light on their possible behavior, especially, at the phenomenologi- 
cal intermediate energy region, where significant non-perturbative effects tend to spoil the 
asymptotic perturbative QCD (pQCD) results like the celebrated quark counting rule that 
predicts the scaling behavior {F, GjxiQ^) ~ 1/Q^ Sl, 13]. Analyses of the pion form factors 



convincingly show that the standard pQCD with only twist-2 effects are much too small 
to explain the currently available experimental data at low- and moderate-energies. This 



calls for the inclusion of non-perturbative corrections from the genuine "soft" QCD [8l-ll2| 
and the sub-leading twists that can give rise to unnaturally large contributions at moderate 
range of Q^-values. in particular, twist-3 enhancements were seen to be quite large in the 
previous studies for the space-like pion form factor 13l-ll7l|. th e sp ace-like kaon form factor 



171 . Il8| . and in the studies of -B — )■ tt transition form factors 19N22|. In fact, the analysis 



presented in J23| shows a possible scenario where the contributions from the twist-3 terms in 
the time-like region can turn out to be exceptionally large. This seemed to resolve the bulk 
of the existing experimental discrepancy between the space- and the time-like pion data. 



In this paper, following 



13, 



23| we extend the analysis to the space- and the time-like kaon 
form factors, where in addition to the twist-2 and twist-3 terms we explicitely include twist-4 
corrections. Thereby, we show that the large twist-3 contributions are indeed a non-trivial 
aspect of our result in comparison with the other twist contributions, e.g., the 2-particle 
twist-4 contributions are explicitly shown to be about a third of the magnitude of the twist- 
2 terms. The paper is organized as follows: the second section briefly reviews the theoretical 
background, the third section deals with the details of our numerical analysis, results and 
discussions of the essential features of our results, and finally, we give our conclusions. For 
the purpose of book keeping, we provide a collection of relevant mathematical formulas in 
the appendix. 

II. HARD AND SOFT KAON FORM FACTORS 

A. Factorized pQCD 

The basic definitions of the space- and time-like electromagnetic form factors are given 
in terms of the following local matrix elements of the electromagnetic quark currents J^™: 

e(P' + P),Fk{Q') = {K^{P')\J-mK^{P)) ; 
eiP'-P),GKiQ') = {K^iP')K-iP)\j;^iO)\0); 

Jir = Y.^i'in,<if. (1) 

/ 



where e is the electronic charge and / is the flavor of the valence quark qj with charge 
e/. In terms of light-cone co-ordinates, P = (Q/a/2, 0, Ot) and P' = (0, Q/\/2, Ot) are the 
incoming and outgoing external kaon momenta in the Breit-frame. In the space-like domain, 
g2 = (p' _ py = _g2 ^ Q^ whereas for the time-like domain g^ = (P' + Py = Q^ > 0. 
Here, Q is assumed to be much larger compared to the kaon mass rriK, so that P and P' 
almost lie along the light-cone directions. 

In our approach, the total contributions to the charged kaon form factors come from 
the factorizable "hard" parts {F,G}^^'^{Q^) calculable within a perturbative framework, 
and the non-factorizable soft parts {F, GY^^^{Q^) relying on some non-perturbative tech- 
niques. The calculation of the hard parts rest on the essential assumption that at suitable 
high energy scales, the form factors are factorizable, i.e., separable into parts dominated 
by short- and long-distance dynamics. The short- distance dynamics are represented by 
the kernel of interactions between highly off-shell partons, above the so-called factorization 
scale ^F- While, the long-distance dynamics below the factorization scale are implicitly 
encoded within the kaonic wavefunctions/distribution amplitudes (DAs) with near-on-shell 
partons. Note that due to the well-known impulse or frozen approximation applicable for 
all high-energy exclusive mechanisms, the dominant contributions come entirely from the 
leading order (LO) Fock state, i.e., a qq valence quark configuration. The higher Fock states 
are neglected with contributions relatively suppressed by higher powers of l/Q^. Fig. [1] 
shows two representative Feynman diagrams (there are 4 diagrams each for the space-and 
time-like cases) with LO hard kernels each having a single hard gluon exchange. These are 
convoluted with the incoming and outgoing kaon DAs to obtain the hard factorized kaon 
form factors. In this analysis, we calculate {F,G}^'^'^{Q'^) up to twist-4 accuracy in the 
two-particle sector, including explicit transverse momentum (/ct) dependence (TMD) of the 
constituent valence partons. The non-factorizable soft contributions, on the other hand. 



can either be calculated using Drell-Yan-West type of wavefunctions overlap ansatz 2J| , or 
from QCD sum rules (QCDSR) incorporating local quark-hadron duality (LD) principle |8|. 
Both these approaches to parameterize the genuine soft contributions are known to give very 
similar results. In this work, we follow the latter approach. The above assertion can then 
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FIG. 1: LO Feynnian diagrams in pQCD for hard contributions to the charged kaon form factors 
in the space-hke (left) and the time-hke (right) region. 

be summarized by 

{F, G}Km = {F, GYfiQ') + {F, G}r^(g2) ; 
{F, Gj^^^'iQ^) = A{F, Cy^'^'^Q^) + A{F, GY^'^'^Q^) + A{F, GY^'^'^Q^) . 



(2) 



The principal inputs for determining the factorized kaon form factors are the coUinear /light- 
cone DAs which encode all the non-perturbative physics. They are universal in nature (frame 
or process independent), in a sense, once they are determined at a certain process, they could 
yield predictions for another. To next-to-leading order (NLO) in conformal twist expansion 
there are one 2-particle twist-2 DA (J)2;K with an axial- vector structure, two 2-particle twist- 
3 DAs - 03.^ with a pseudo-scalar structure and ^g.^ with a pseudo-tensor structure, and 
finally, two 2-particle twist-4 DAs A4-k and M^-k = gi;K — (p2;K both having pseudo-scalar 



structures 
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-|27| . As an example for K , we display the twist-2 DA in terms of the following 



pseudo-scalar matrix element with ^ = 2x — 1: 

{0\u{z)-f,-f,s{-z)\K-{P))=tP, [ cixe^«^^)02;i^(a;,/x') 

with the normalization condition 

Ik 



N, 



2:K 



^-'-^^^^'^'^- 2V2N: 



where fx is the kaon decay constant defined in the local limit 2; — ?■ by 

{Q\u{Q)j,^,siO)\K-{P))=tfKP,. 
5 



(3) 



(4) 



(5) 



In the above equations, x is the colhnear/hght-cone momentum fraction (xj = k~^/P~^) 
carried by the individual valence quarks {x for the s quark and x = 1 — x for the anti-quark 
u.) Note that the gauge-connection factor in the above matrix element is assumed implicitly. 
To the leading logarithmic accuracy 02;X satisfies the well-known ER-BL evolution equation 
p, l6| and can be expressed as an irreducible representation of the special collinear conformal 
group SL(2,R), in terms of standard Gagenbauer Polynomials Cj (0- 

<P,M^,f.') = 4S(x) X:«n'(/^o)C^r(0 (^J + 0{a,) , (6) 

with the asymptotic twist-2 DA given by 



,(as)/ \ J / 2 > \ "^JK 



X) = 4>2;k{X-, /i^ — )■ Oo) = X(l — X) . (7) 



The standard QCD MS running coupling as(/i^) to two-loop accuracy is given by 

«.(/i^)_ 1 /3iln(ln(/xVA^CD)) 



(8) 



/3oln(/xVA^CD) /3o'ln'(^VA|cD) 
with Aqcd ~ 0.2 GeV, (3q = (llA'c - 2A^/)/12 = 9/4 and /3i = (51 A'^ - 19A^/)/24 = 4 for 
Nc = Nf = 3. The ratio of the QCD couplings represents the renormalization group (RG) 
evolution of the Gagenbauer moments a^ from the normalization scale /iq ~ 1 GeV to the 
generic scale fi, with LO anomalous dimensions given by 



4 fl ^1 



The Gagenbauer moments represent the genuine non-perturbative inputs to the DAs and 
are usually determined using lattice simulations or from light-cone sum rules (LCSR). In 
this work, we use the latter inputs, since the moments for the higher twist DAs are yet to be 
determined precisely in Lattice QCD. Note that the lower order moments in both approaches 
are known to be in good agreement with each other. However, dealing with such an infinite 
number of terms in the non-asymptotic DAs become a matter of technical challenge as the 
higher order moments are extremely difficult to determine. Hence, for practical simplicity of 
calculation, one truncates the DA series up to the first couple of terms only. Moreover, the 
increasing anomalous dimensions tend to suppress the higher order terms. In this analysis, 
we consider the series up to the term with the second moment af . The rest of the non- 
asymptotic collinear DAs, i.e., the 2-particle twist-3 and twist-4 DAs which we also consider 



in this work, have more elaborate expressions and are, therefore, relegated to the appendix 
along with their RG evolutions. A summary of the relevant DA parameters determined from 
LCSR at the normalization scale of /io ~ 1 GeV is presented in Table. [H 

A common feature of DAs derived from LCSR is that they are endpoint dominated 
due to large kinematic enhancements when the light-cone momentum fractions tend to the 
endpoints (i.e., x — )■ 0, 1). One possible way to suppress such artificial enhancement is to use 
the Brodsky- Huang-Lepage (BLH) Gaussian parameterization 28|, where the intrinsic TMD 
of the valance partons within the full kaon wavefunction "^t-^K (for each twist t = 2,3,4) is 
explicitly modeled by including an addition wavefunction S^-x, i.e.. 



^f;i^(x, kr, f/, M{u,d,s}) = A;K 0t;x(a;, /i^) ^^.^(x, k^, M{u,d,s}) , 



(10) 



where the form of I^t;K is chosen similar to that of a harmonic oscillator wavefunction that 
can maximally suppress such endpoint effects and given by 



Et.K{x,kT,M{u,d,s}) 



16vr^A^;X 


o2 


(Ml + ]^ , Ml, + \.l\-\ 


a;(l — xj 


-K;K 


-^ — - + ^ 

\^ X 1 — X J 



:ii) 



Note that the constituent quark masses M.{u,d,s} are introduced to parameterize the QCD 
vacuum effects, while the parameters A^k and (5t;K for the individual twists are phenomeno- 
logically extracted as described in the next section (also see, 17|). Next to obtain the 
full TMD modified kaon DAs Vf^K in the impact parameter or 6-representation, we use the 

5[, yielding 

1 «, 



Brodsky-Lepage definition of the DA 



IGtt^ 



*t;M(a;,kT,Ai^A^{„,d,4) 



At-K4>t;K{x,^i^)ey:.V 



-A 



t;K 



A^2 M 



uA 



X 



1 — X 



xexp 



h'^xil 



X] 



4A' 



t;K 



(12) 



In Figs. [2] and [3l we show the various collinear DAs (which are endpoint enhanced) and 
the modified BHL DAs (which are endpoint suppressed) respectively. We also display the 
corresponding asymptotic forms of the DAs. 

The inclusion of the TMD in the hard scattering kernel at the same time also serves 
as a natural regulator for possible endpoint enhancements. However, this leads to the 
appearance of large logarithms in the kernel due to incomplete cancellation between soft 
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FIG. 2: The collinear twist-2 and the 2-particle twist-3 and twist-4 DAs for the kaon (modulo, 
the normahzations Nt-x), shown using sohd (red) hnes, i.e., (a) 02;i<'(a^5/^o)' (b) 'Ps-xi^^ f^o)^ (^) 
4>^.j^{x, fj,Q), (d) mj^Ai-xix, Hq), and (e) rnj^gi-xix, fJ-Q), along with their respective asymptotic 
DAs, shown using dashed (green) lines. The DAs are defined at the scale fiQ = 1 GeV. 

gluon bremsstrahlung and radiative corrections that may spoil the perturbative convergence 
and, hence, the validity of the collinear factorization. While, the large single-logarithms 
such as as In Q^ can be effectively tackled using usual RG techniques like UV divergences, 
the large double-logarithms or Sudakov logarithms involving TMD such as a^ ln^(Q^/A;|,), 
arising from the overlap of the leading soft and collinear kinematic regions of radiative 
gluons, can not be similarly handled in ordinary fixed order perturbation theory. The 
alternative is to use resummation techniques to all orders in the strong coupling constant as 
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FIG. 3: The BHL modified twist-2 and the 2-particle twist-3 and twist-4 DAs for the kaon (modulo, 
the normalizations Ni-^k)-, shown using solid (red) lines, i.e., (a) p2;Kix, fJ^^), (b) P^.j^{x, fi^), (c) 
■P3'!'^(x,/ig), (d) m'j^Vfj^{x,fj,Q), and (e) 'm'j^V^.j^{x, ij.q), along with their respective asymptotic 
DAs, shown using dashed (green) lines. The DAs are defined at the scale /ig = 1 GeV. 



which organizes the double-logarithms within exponential Sudakov factors to get eventually 
systematically absorbed by a re-definition of the DAs. Such Sudakov factors represent the 
perturbative tail of the DAs and suppress non-perturbative enhancement that arise from 
constituent partonic configurations which involve large impact space separations. For a 
review of the Sudakov form factors and their application to exclusive physics, the reader is 
referred to 
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m 



29 



30| . There may be other radiative coUinear double- logarithms such as 



tts In X which may be resummed using threshold resummation 20|, |2l| to suppress additional 



collinear enhancements in the kernel. The threshold resummation along with the Sudakov 
resummation arising from different subprocesses in pQCD factorization provides natural 
suppression to the endpoint and other non-perturbative enhancements and are relevant in 
the range of currently probed energy /momentum transfer values. The upshot is that the 
hard perturbative contributions are enhanced relative to the non-perturbative contributions 
improving convergence significantly and making pQCD evaluation of exclusive form factors 
self-consistent toward lower values of Q^, where it may not be otherwise justified. 

Such techniques of systematic organization of the potentially large logarithmic contribu- 
tions is a modification from the standard collinear factorization and is generally termed as 
the " kf-factorization" , that has been widely applied to inclusive as well as exclusive pro- 



cesses 3l|]. However, unlike the familiar collinear factorization theorem, the fc^- factorization 
is currently considered only at the level of a conjecture which is yet to be proven to all orders 
in perturbation theory (this is a highly debatable issue and, in fact, not yet fully recognized, 
e.g., see [32] for a different viewpoint.) To demonstrate that /cT-factorization is indeed a 
systematic tool, demands higher order calculations which may be very challenging. How- 
ever, in this paper, we shall implicitly assume the validity of such a modified factorization 
without proving it and restrict ourselves at the tree level analysis of the /cT-dependent hard 

the fcy-factorization was proven at the level of twist-2 accuracy. 



kernel. Moreover, in 



while the collinear factorization was explicitely shown to be valid at the twist-3 accuracy 
in the case of the 717* — )■ 7 transition form factor. Our analysis is, therefore, based on the 
key assumption that the same formalism could be straightforwardly extended to the elastic 
kaon form factors. 

At the leading order ~ 1/Q^, the twist-2 and the two-particle twist-4 terms contribute 
to the hard kernels which have exactly the same expression given by 

^(t=2,4;LO). ^21 , 2n _ ±16n Cf asiii^) xQ"^ 

\ard (X,y,g,k^,lT,/i)-^^^2±k2^)(^^g2±(k^_l^)2)' (l^j 

while, the C(l/Q^) power suppressed 2-particle twist-3 and twist-4 hard kernels are, respec- 
tively, given by 

3271 Cp as{fi'^) X 



^hard (^5 1/5 <5 5 ky. It, /i ) /„^2 _L 1,2 N /_. ^2 _L /U 1 ^2^ 



(xg2±k2,)(a;?/g2±(kT-lT) 

rpl t=4;lM^UJ / ^2 1 1 2\ 4:OTT Lp (Mgyfl j ft a\ 



^(t=4;NLO)/ ^2 1 1 2^_ jSTtXpOjjJ^'^ 

^ hard la;, y, t^ , Kt, It, /i j - ^^^^ ^ ^^^^ ^^^ ^^ ^ ^^^ _ ^^^ 
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In the above expressions, Cp = 4/3, the "+" signs correspond to the space-hke case and 
the " — " signs correspond to the time-hke case; k^ and 1^ are, respectively, the initial and 
final relative transverse momenta of the valence quarks, and x and y are the corresponding 
light-cone momentum fractions. Note that the factors in the denominators that arise from 
the parton propagators develop poles in the time-like region. 

To obtain the hard form factors, we use the following momentum space projection oper- 
ator for the DAs with the different twist structures: 

iK _ M itL. /,T, ^_2 ,T,A a2 \ , _2 . / f ^ 



- f^Kl, l^n,K - ^^M.^'^^'^^^3;X + l^.^P" n,K ^fc^) } ', P = \P\n , (15) 

where hk = ^ J^^ GeV is the " chiral-enhancement" parameter arising in the standard 
definition of the 2-particle twist-3 DAs (see, appendix), "^^.^ = '^I.j^ — '^2;K, (^Ut — 'j'^d/dkl^, 
n = (1,0, Ot) the unit vector in the "+" direction, and n = (0,1,0^) the unit vector in 
the "— " direction. Setting the renormalization/factorization scale to the magnitude of the 
incoming or outgoing kaon momentum i.e., jj, = \P\ = \P'\ = Q/\/2 and convolving the 
projection operators for the kaon DAs with the hard kernels using factorization formula 
(symbolically, A^K;out ® ^im?d ® Mk-m), we have 

(P'±P)^{F,G}r'(Q') 

= Jo ' y T^ T^ [ — 3 — ) •^"p ^-'^^ ■ ^' - '^^ ■ ^^] 



Tr 



YMk-ouI^i^ knternan^, MK;in ^ ,. 

■ + 3 diagrams 



WRGE(t,/i) 



X St{x) exp [-S{x, y, |kr| ~ l/&i, |1t| ~ 1/^2, /i)] , (16) 

where kintemai and kg are the internal quark and gluon momenta, respectively, as shown in 
Fig.l. Also, in the above equation 



WRGE(t,/i) = exp 



4 / ^7g(as(;u2)) 
t /" 



;7>.(^^)) = -^^^, (17) 

71 



represents the RG evolution factor for the scattering kernel from the "upper-factorization" 
scale t = raa,x(y/xQ, l/6i, I/62) to the renormalization scale /i = Q/^/2, and 7^ is the quark 
anomalous dimension. The expression for the Sudakov exponent S{x,y,bi,b2,Q) (after 



11 



absorbing the RG factor from the kernel) is given by [3 

S{x, y, bi, 62, Q) = s{xQ, 61) + s((l - x)Q, 61) + s{yQ, 62) + s((l - y)Q, 62) 
+2 f ^7.(«s(/^^)) + 2 r ^7>.(/i^)), 



where 



s{xQ, 1/6) = s(x/i, 1/6) 



XfJ, 



d^ 



i/b ^^ 



\n\ ^\ A{as{ji'')) + B{as{ji')) 



(18) 



(19) 



where the "lower- factorization" scales I/61, I/62 > Aqcd serve to separate the perturba- 
tive from the non-perturbative transverse distances which are also typically the scales that 
provide a natural starting point of the evolution of the kaon wavefunctions. In the above 
equations, the so-called "cusp" anomalous dimensions A and B, to one-loop accuracy are 
given by 



A{as{^l^)) = Cf 



asifi^ 



+ 



TX 



67 TT^ 
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BM.=)) = f^ln(^ 



27 9 y 27 -^ 3^ V 2 

27E-1 



3 7B 



as{ii^ 



TX 



(20) 



The exact form of the threshold resummation "jet" function St[x) in Eq. [T6] involves a one 
parameter integration, but in practice it is more convenient to take the simple parameteri- 



zation proposed in 



mm- 



St(x) 



2i+2T(3/2 + c) 



[x{l-x)Y , 



(21) 



0Fr(i + c 

where the parameter c ~ 0.3 for light pseudo-scalar mesons like the pion and the kaon. 

Now we present the factorized result for the hard kaon form factors up to twist-4 correc- 
tions as follows: 

{F, G}5,^'-'^(g2) = S{F, G}f'''\Q^) + 6{F, G}f''''{Q') + S{F, G}™(g2) ; 

S{F, GYr'^Q') = HF, G}*/'^*'''^°(g') + S{F, G'}*/^*"'''^°(Q2) , (22) 

where the leading twist-2 and twist-4 corrections are expressed by the following integral 
representations in the impact parameter space: 

/■I /-oo 

<5{F, G}*/^*2 (g2^ ^ ^1^^ ^|twi.t4;Lo^g2^) _ s2ttQ^Cf / X dxdy / 61^6162^62 as{t) 

Jo Jo 

- 1 52 

X ±- V2;k{^, 61) V2;K{y, 62) Tmj,^ V2;k{^, 61) vt^Ay, h) 

T m^ I Pt^x, 61) P^-MV^ h) + C^ «&i, ^^t) 

o 

X i:f±(x, y, Q, 61, 62) exp [-^(x, y, 61, 62, Q)] 'S't(x) , (23) 
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while, the power suppressed twist-3 and twist-4 corrections are given by 



6{F, G'}*/^*3(Q2) + 6{F, G}™'^^°(Q2) = 32nQ^CF / dxdy / hdhb2db2 a,(t) 



/^i 



^[xVi.^A^M)Vi.^^{yM, 



'^^""^ i-KK{^^^^)nAyM) 



'-,K,Ki^^b,)Vl^^iy,h) + 



3m 



K 



2Q2 



6 dx 



dCVlACM 



X H±{x, y, Q, bi, 62) exp[-S'(x, y, 61, 62, Q)] St{x) , 



(24) 



where Vt;K{x,b) = Vt,K{x,b,l/b'^,M{u,d,s}), and VIk = VIk - ^2;X- Note that the "LO" 
and "NLO" used in the above equation should not be confused with the usual terminologies 
associated with perturbative expansions in terms of as but rather in the sense of operator 
product expansion (OPE) terms. In the impact representation, the space- and time-like 
hard kernels (the part of the scattering kernel that is common to all the twists) could be 
expressed in terms of standard Bessel functions Kq, Iq, Jq and Hq and are given by 



H+{x,y,QMM) 





KQ{^/xyQb2) 

X [0(61 -62)i^o(v^Q&i)/o(v^Q&2) 
+ e{b2 - 6i)Ko(v^g62)/o(v^Q&i)] 



(25) 



H_{x,y,Q,bi,b2 



f ) i^i"' 



xyQb2) 



e{bi - b2)H^'\V^Qbi)JoiV^Qb2 



+ d{b2 - bi)H^o\V^Qb2)MV^Qbi) 



(26) 



where if+ is a real- valued function and H_ is a complex-valued function of real arguments. 
Apropos of our derived formulas Eqs. |23]and[2U it is noteworthy to mention that in 34 1 
it was suggested that the Sudakov factors must be analytically continued from the space- 
like to the time-like case. This may not be generally true. The Sudakov factors in [35 1 
(see, section 3.1 of this reference) arise directly from "form factor- type" kernels, which are 
not universal quantities and may vary with processes. There the analytic continuation is 
perfectly justified. However, for an approach based on the factorization theorem, one uses 
"universal" Sudakov factors S{Q) arising from the overlap of the soft and coUinear processes 
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below the factorization scale, as in the present context. As explained in [36|, these Sudakov 
factors are to be considered as an integral part of the DAs and, thus, they are universal 
quantities as well, depending only on the magnitude of the energy scale Q > 0. Note that 
the Q dependence of the Sudakov factor in Eqs. [TH]and[Tni stems from the dependence on the 
collinear components of the external pion 4-momenta which are given by P+ = P'^ = Q/\/2, 
in the Breit-frame. As such, it is important that one does not analytically continue but rather 
use the same Sudakov factor in both the space and time-like cases. 

B. Non-factorizable Soft QCD 



In 



a 



it was shown that the space-like low-energy pion data below Q^ ~ 10 GeV^ is 
dominated by the soft pion form factor which accounts for more that 70% of the data. 
Such soft QCD contributions are non-factorizable and are beyond the realm of ordinary 
perturbation theory. Since, no systematic method is currently available to calculate these 
non-perturbative effects, one is compelled to use some model ansatz to obtain a rough 
estimate of their contributions, viz, in J8| the soft pion form factor in the space-like region 
was calculated using the Local Duality (LD) model in QCDSR. In our present work, we 
extend the same result to the space-like kaon form factor which is then given by 

Fk (Q)lLD = l-(i + 4,^/g2)3/2-gI + ^(g^)'^o-4vr/^. (27) 

Now, on one hand, an ah initio derivation of the corresponding time-like soft form factor 
seems a priori unfeasible using QCDSR, since the LD principle is strictly applicable for the 
space-like region only. On the other hand, a naive analytic continuation of the space-like 
formula, i.e., by a replacement of Q^ — ;■ —Q^-, leads to an undesirable pole in the denominator 
of the soft form factor: 

^soft ,q2. 1 _ 1 - 6^o/^' ^ ^ , oi—) r28l 

^i.;analyticW )^^ (^ _ 4,^/g2)3/2 ~ g4 + ^^ qJ ' i^«J 

Since, here our primary goal is to obtain an estimate for the smooth continuum part of 
the kaon spectra for intermediate energies which in reality is, however, dominated by low- 
energy time-like resonances that obscure the smooth continuum. With the "over-simplified" 
assumption that these resonance peaks behave as background "noise", superimposed on a 
smooth continuum spectrum, we choose the functional form of the time-like soft form factor 
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to be the same as that of the space-hke expression, which is a smooth function for the entire 
range of Q^ we consider, i.e., 

grW-)^l- (,Vl':;^^;. .0(^). (29) 

Moreover, for large enough Q^ ~ above 5 GeV^, both expressions Eqs. [27] and |2S] when 
expanded in inverse powers of Q^ yield the same leading term of (9(4i). Hence, the particular 
choice of the soft form factors should not matter significantly at large-Q^ values where the 
perturbative predictions become more reliable and dominant. 

In the present context, a vital aspect deserves some consideration. Since, the inclusion 
of the soft form factors has been somewhat ad hoc^ without any correspondence among 
the hard and the soft contributions, there could be chances of possible double- counting of 
contributions especially at low energies. Thus, it becomes clear that we must correct the 
hard factorized results in the low-Q^ region to ensure that the respective contributions lie 
within their domains of validity. This is achieved by enforcing the gauge invariance condition 
through the vector Ward-identity {F, G}k{Q'^ = 0) = 1, which is a priori not ensured in 
perturbative calculations. Since the soft form factors satisfies {F,GYk^{Q^ = 0) = 1, we 
must have {F, Gy^'^'^{Q^ = 0) = 0. But this is unfortunately not satisfied by Eqs. [23] and 
[21] where the contributions tend to diverge rapidly in the vicinity of Q^ = 0. Therefore, the 
essential task is to match the large-Q^ results of {F,G}^^'^{Q'^) with the low-Q^ results of 
{F,GY^^^{Q^). Here we shall modify the argument given in pjj] for the twist-2 case to be 
applicable for the twist-3 and twist-4 power corrections. The simplest way is to "power- 
correct" for the singular ~ l/<5^ (leading twist-2 and twist-4) and ~ ^/Q'^ (sub- leading 
twist-3 and twist-4) behaviors, respectively, at small Q, by introducing some characteristic 
low-energy mass scale Mq that may lead to the onset of the genuine non-perturbative soft 
dynamics. For the soft form factors modeled via LD principle, the scale Mq = 2sq is a 



natural choice 37(]. It can then be shown that for the leading twist-2 and twist-4 hard 



corrections, it is sufficient to make the modification [12| 

6{F, GY^'^'^Q^) + 6{F, G}*/^*''^°(g') 

^ A{F, G}™(Q2) + A{F, G'}*/'*''^°(Q') 

= (^^^^) (^{^' GYr'^iQ') + HF, GYr'^'^^'^iQ')) . (30) 

For the case of the sub-leading twist-3 and twist-4 hard corrections, we perform the following 
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replacement: 

,, o .N - — ■ (t=3,4) /Lf4 , — . ■ {t=3,4) /\/f4 

5{F, G}r '') m = 5{F, G}^. (g^) ^ ^ 5{F, G}^ (g^) ^^^^ , (31) 

where we write (5{F, G}i^=^'^^ = 5{F, G}^'^'3^(5{F, G}^^'^^^'^^^ for brevity. Now, to maintain 
the Ward-identity, we correct for the wrong g^ = hmit of the above expression 

/^ o .^ ' (t=3,4) -—- {t=3,4) /Vf4 

A{F,G}J=^'^Hg2) = -5{F,G}^ (g^)$n(gVMo^) + <5{F,GK. ^^'^M^Tg^' (^2) 

where we introduce the smooth function $„(-2) (with z = Q'^/Mq) with the essential property 
that $„(0) = 1 and z^n{z) — ;■ as ^ — )■ oo, for a suitable choice of the positive integer 
n, to preserve the asymptotics of S{F,G}j^~ ' (g^). A natural choice for $n(^) could be 
$„(z) = 1/(1 + z")^, concurrent with the ~ 1/g^" scaling behavior of the respective power 
suppressed terms. For the present purpose, it is sufficient to take n = 2, yielding 

A{F, G}^ (g ) - S{F, G}^ (g ) ^^^^-^^ \^1 - j^p^,) 

= S{F,G}^t'^'\Q^)(-^)\ (33) 



M4 + g4 

In principle, this can also be achieved with larger integer values of n that would lead to 
(g^"/(^o" + Q^")) ill front of the hard parts. However, as n — )■ oo, this factor becomes a 
step function, which is no longer smooth. Hence, a minimal value of n is preferable and we 
arrive at the Ward-identity modified result: 

s{F, GYr'HQ^) + HF", G}*7^*''''^°(g2) 

^ A{F, GYr'HQ^) + A{F, G}*/^*^'''^°(g2) 

= (i^l^) {HF, GYr'HQ') + 6{F, G}r*^^^^°(g2)) . (34) 

The pre-factors only alter the low-energy behavior of the hard contributions and ensure 
the correct power-law in maintaining a smooth matching between the large-g^ behavior of 
{F, G}^j^"^{Q^) to the low-g^ behavior of {F, GYg\Q^) (see. Fig. H). This leads to our final 
expression for the total electromagnetic kaon form factors, correct up to C(^) accuracy, 
given by 

{F, G}k{Q^) = {F, GY°'\Q^) + A{F, G}™(g2) ^ ^^p^ ^|twist3^g2) ^ ^1^^ (<.|twist4(g2) ^ 
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where 



{F,G}r{Q') = l 1 + 6WQ^ 






A{F, G}r*'(g^) = ( ^-^ ) 6{F, GYr'^iQ 

twists /,n2\ _ / ^ \ xr I? /^l twists /,n2\ 



A{F, G}™(Q2) = ( ^^^^) H^, GYr'''''^''{Q') 

^ \ i: r r:i ^t twist4;NLO /^2 



5{F,G'}™'^'^^(Q"), (35) 



where the (5's are replaced by the A's to include the respective pre-factors. 

III. RESULTS AND DISCUSSION 

To obtain the BHL Gaussian parameters of the kaon wavefunctions, we use the following 
two sets of constraints valid for the individual twists (t = 2, 3, 4): the first set of constraints 
is obtained from the leptonic decay K ^ fi + u^^, and given by 

dx / — -^ '^t;K{x, kr, fll, M{u,d,s}) = Nt;K , (36) 

with Nf-K being the normalization constant for the collinear DAs; and the second follows from 
the phenomenological fact that the average transverse momentum of the valence partons in 
light mesons is about {kT)jKn--- ~ ^.35 GeV, i.e., 

9 

Jdx J d'^kr I ^t;x(a;, ky, /xg, M{u,d,s}) \ 

The Gaussian parameters determined in this way for /xq ~ 1 GeV are collected in Table. [TTl 
Note that due to the rather mild scale dependences of these parameters, which practically 
remain constant for the entire range of intermediate energies that is considered in this work, 
their scale variations have been kept fixed to reduce the numerical complexity. However, we 
do consider their variation with the changes in the collinear DA parameters, summarized 
in Table. HI that is required for our estimation of the theoretical error. Once all the phe- 
nomenological parameters are determined, we proceed to calculate the hard contributions. 
For calculations, we use the full non-asymptotic collinear DAs derived from LCSR 25|-|27|. 
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K^ parameters 


At ^0 = 1 GeV 


units 


K^ parameters 


At ^0 = 1 GeV 


units 


mu,d 


5.6±1.6[27] 


MeV 


4 


0.25±0.15[27] 


- 


rus 


137 ± 27[27] 


MeV 


Ik 


1.22/^, U = 131[381 


MeV 


Mu,d 


0.33 


GeV 


fsK 


0.0045 ± 0.0015[27] 


GeV2 


Ms 


0.45 


GeV 


^3K 


-1.2±0.7[27j 


- 


rriK 


493 


MeV 


4 


0.20 ± 0.06[27j 


GeV2 


«f 


0.06 ± 0.03[27] 


- 


^4.K 


0.2±0.1[27] 


- 



TABLE I: Various input parameters for twist-2, twist-3 and twist-4 wavefunctions. 



At,K{t = 2,3,4) 


At ^0 = 1 GeV 


units 


(/3i;^)2(t = 2,3,4) 


At /io = 1 GeV 


units 


M;K 


2.06(2.07)as 


- 


{h,K? 


0.78 (0.89)as 


GeV-2 


A'P 


2.23(2.28)as 


- 


m,K? 


0.71 (0.79)as 


GeV-2 


^3;K 


2.08(2.07)as 


- 


^K? 


0.88 (0.89)as 


GeV-2 


aA 


2.07(2.00)as 


- 


Kk? 


0.91 (0.93)as 


GeV-2 


A^ 


5.22(2.28)as 


- 


^K? 


0.57(0.79)as 


GeV-2 



TABLE IL The phenomenologically determined Gaussian parameters for twist-2, twist-3 and twist- 
4 wavefunctions. The numbers in the parentheses (• • • )as correspond to parameters for the asymp- 
totic wavefunctions. 

In the Fig. H we plot the individual terms of Eq. |35l i.e., {F, G}"/S A{F, G}*/'^*^ 
A{F, GYit^^'^ and A{F, G}^'*^*^, which should give an idea about the relative magnitude 
of each contribution for intermediate values of Q^ up to 30 GeV^. For comparison, we also 
display the results obtained without including the pre-factor modifications, which do not 
show any appreciable difference for Q^ values beyond ~ 5 — 10 GeV^. As expected, the 
standard twist-2 contributions are much smaller compared to the soft QCD and the twist-3 
power corrections at moderate-energies. However, the twist-4 contributions are seen to be 
indeed small (about 1/3 of the magnitude of the twist-2), which are, in fact, negative in 
the space-like region. In the time-like region, since all the hard contributions are complex, 
it only makes sense to plot the modulus of the individual twist corrections. It is notable 
that the general enhancement of all the time-like hard contributions relative to the space- 
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e'8[f,|G|];"'^'2 
e'A[f,|G|];"'^" 

q' [f ,|g|]/°" 




FIG. 4: Relative contributions of the soft {F,G}^^* (double-dot black lines), twist-2 A{F, GY/'''^'^ 
(thick double-dot dashed orange lines), twist-3 A{F, G}^"^*'^ (thick dotted blue lines), and twist- 
4 A{F, GY^^^^^ (thick solid red lines) terms in Eq. [35j The same terms without the pre-factor 
modifications are also displayed. 



like ones can be attributed to the time-like parton propagators developing poles that are 
absent in the space-like region. To illustrate this point, it is useful to plot the part of the 
hard kernel H±{x,y,Q,bi,b2) that is common to all the twist corrections to the hard form 
factors. Fig. [5] shows the variation of the space- and time-like kernels H± (in the impact 
parameter space) as a function of Q^ for some arbitrary fixed values of the parameters 
x,y,bi and 62- It immediately becomes clear that the real- valued space-like kernel H^ has 
a rapidly decaying exponential behavior, whereas the complex-valued time-like kernel if_ 
has rather large amplitude oscillatory real and imaginary components which decay very 
gradually with increasing Q^. In reference to Eqs. [T3]and[T4l we note that if {x, y} ^ 1 ^ 
{xQ"^, xyQ"^} ~ {ky. If.} < Q"^, and if {x,y} ~ 1 ^ {xQ"^, xyQ"^} ~ Q"^, so that the terms 
in the denominators tend to cancel each other in the time-like but not in the space-like 
domain. This explains why the amplitude of the time-like oscillations in if _ grow larger and 
larger near the endpoints x, y — )■ 0. 

The most striking feature of our results in Fig. H] is the anomalously large twist-3 contri- 
bution in the time-like region, similar to what was seen for the pion 23| , dominating all the 
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0.5 



-0.5 




H^{x, y. or. b„ b^) 

\H_(x, y, d, b,, b2)\ 

\m[H_(x, y, q\ b„ b2J\ 

Re[H_(x, y, ^, bj, b^)] 



X = y = 0.1 

bi = h = 0.6/Aecfl 



10 



20 

Q^ [GeV^] 



30 



40 



50 



0.5 



0.25 



-0.25 



-0.5 



H^{x, V, e , b„ b2) 

\H_(x, y, (i, b,, b2)\ 

\m[H_(x, y, q\ bp b2)] 

Re[HJx,y,Q\bj,b2)] 



M. 



x = y = 0.9 



h, = b2 = 0.6/Kqc^ 



10 



20 



30 



40 



q2 [GeV^] 



50 



FIG. 5: The space- and time-like hard kernels Hj-^x, y, Q, bi, 62) in the impact space representation 
for two sets of choices for the collinear momentum fractions with arbitrary fixed 61, 62: x = y = 0.1 
(left plot) and x = y = 0.9 (right plot). 

other corrections for the entire range of low- and moderate-energies. This huge asymmetry 
between the space- and time-like twist-3 contributions comes from the additional parametric 
enhancement of the twist-3 DAs due to the chiral parameter fiK which makes them partic- 
ularly sensitive to the chiral scale. It is the combination of this parametric enhancement 
along with the occurrence of the time-like poles in the hard kernel that leads to such a 
characteristic anomalous twist-3 behavior which is completely missing in the twist-2 or even 
in the twist-4. At this point, one may also worry about possible large contributions from the 
3-particle twist-3 sector (related to the 2-particle twist-3 sector through QCD equations of 
motion) that was not considered in this work. Here, we note that such a possibility can safely 
be precluded since the 3-particle twist-3 DA receives large parametric suppression from the 
non-perturbative parameter f-^x ~ 0.0045 GeV^, numerically very much smaller compared 
to the analogous 2-particle twist-3 parameter fix ~ 1-5 GeV which greatly enhances the 
contribution from the 2-particle sector. 

Further, it is important to note that (a) the "active" soft gluons that may also arise 
from the 3-particle twist-3 DA or higher twist DAs likewise, bring about additional power 
corrections and, therefore, can be safely neglected at large-Q^ values, and (b) the "long 
distance" soft gluons that may be a possible source of the breakdown of TMD factorization, 
can not probe the small "color-dipole" configurations of the qq hadronic bound state at 
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q^ [GeV^] 



FIG. 6: The total scaled kaon form factors, denoted by the thick solid (red) lines; the soft form 

Q 

factors {F,GYk^^ denoted by the double-dot (black) lines; and the asymptotic QCD result [7|], 
denoted by the short-dashed (orange) lines. The shaded area is our estimated theoretical error. The 

m n 

experimental data taken from [2, l3|] and the phenomenological result [39|] is shown for comparison. 



high enough Q^ (color-transparency). The remaining collinear gluons are assumed to be 
effectively tackled within the present TMD factorization scheme, where the inclusion of the 
2-particIe twist-3 corrections indeed turn out to be the most crucial aspect at the moderate- 
q2 rggij^e. Note, however, that all such non-perturbative power corrections including the 
soft contributions rapidly fall-off with increasing Q^, and beyond ~ 50 — 100 GeV^ the 
standard twist-2 terms start dominating the asymptotic regime, yielding back numerically 
the bona fide asymptotic behavior given by the Farrar and Jackson result |7|. 



{F,G}^^Q 



^^y ('/O^'i 



ircasiQ')!' 



K 



Q' 



(38) 



Our final prediction for the total scaled electromagnetic kaon form factors {F, G}k (from 
Eq. 1351) up to twist-4 accuracy in the range of intermediate energies/momentum transfers 
is presented in Fig. El along with the result for the soft form factors {F, G}^^* and the 
standard asymptotic QCD result of Farrar and Jackson [7| for comparison. To estimate the 
theoretical error we studied the variation of the wavefunction parameters provided in the 
Tables [T] and [TTl In addition, we varied the chiral parameter fj,K = Tnj^/{mu + mg) which 
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is often taken to be slightly lower ~ 1.3 — 1.5 GeV in the literature |17l . ll9H23l. l26l . |40| . |41 | 
than its naive value about 1.7 GeV expressed in terms of the current quark masses. In this 
analysis, we take fix = 1-5 ± 0.2 and include its variation in the error estimate. The shaded 
area, thus obtained, can be regarded as our rough estimate for the theoretical error, where 
the solid (red) curve corresponds to the central values of the parameters. While our result 
is relatively insensitive to the choice of the parameters in the space-like region, the time-like 
result turns out to be very sensitive to the choice of fix whose variation alone amounts 
for more than 90% of the error-bar. The error due to the rest of the model parameters 
is generously overestimated to include possible uncertainties due to the soft parts which 
we do not a priori take into account. Thus, we should stress that our pQCD based error 
estimate in the low-Q^ region (which apparently looks small) must be considered in a very 
conservative sense and can not be taken seriously below ~ 5 GeV^. A more rigorous error 
analysis is impossible at the moment due to poor quality of the experimental data. 
Several comments are now in order: 

• The width of our error-bar is large enough to completely subsume effects due to further 
inclusion of higher-twists (e.g., twist-5 and twist-6), sub-leading Fock states and higher 
helicity components whose contributions should be tiny, not exceeding even 1%. 

• Our LO (in a^) scattering kernels are apparently gauge dependent arising from the con- 
tributions of the single hard gluon propagator. However, in the context of the 717* — )■ 7 
transition form factor, it can be shown through a systematic order by order calculation 
using AjT-factorization that there is indeed a cancellation of the gauge dependences be- 
tween the quark-level diagrams of the hard kernel and the effective diagrams of the 



pion wavefunction J42|, so that the net result turns out to be gauge-invariant to all 
orders. It is, thus, believed that the same technique can be straightforwardly extended 
to other hadronic elastic and transition form factors, including the present context of 
the kaon form factors, at least up to the level of NLO corrections. 

The factorized hard form factors further suffer from renormalization/factorization scale 
dependent ambiguities that typically emerge from the truncation of the perturbative 
series and would be absent if we were able to obtain an all-order result in the QCD 
coupling Og. To minimize the scale dependence in our present investigation, we adhere 

mk^, as 



to a fixed prescription with the scales set to the momentum transfer Q 
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mentioned previously in the context of the Sudakov factor. In this way, we hope to 
improve the rehabihty and self-consistency of the perturbative prediction and reduce 
the influence from higher-order corrections. 

Nevertheless, a naive estimation of the NLO twist-2 contributions to the kaon form 
factors, using available NLO calculations for the pion form factor in asymptotic QCD, 



can be roughly expressed as J43| 



Q'{F, G}^^" ^ (0.903 GeV^) al{Q') ^ , (39) 



n 



which yields a rather nominal contribution ~ 20 — 30% that is roughly of the same 
order of magnitude as the twist-4 contributions obtained in our analysis. It is to be 
noted that the above estimation is based on the usual collinear factorization approach 
|43l | which did not take TMD into account. Including the k^ dependence of the kernel 
might further reduce the NLO corrections, as was shown in the cases of the pion 



44| and the nucleon [45!] form factors. It goes without saying that a full systematic 
NLO calculation (including twists-3) within the TMD factorization scheme, which is 
missing until now, would be indispensable in resolving this issue about the definitive 
magnitude of the sub-leading corrections. 

On the experimental side, as seen in Fig. Ej currently the space-like region is completely 
devoid of data points at Q^ values higher than ~ 0.2 GeV^. This makes it difficult, if 
not impossible, to compare such low-energy data with our predictions based on a pQCD 
approach which becomes unreliable and diverges rapidly in the vicinity of the Landau pole 
Aqcd ~ 0.2 GeV. For the time-like region, there existed some older kaon data at relatively 
higher energies but with very poor statistics [2j. For such measurements the data above 
Q^ > 4.7 GeV^ had either upper limits or errors > 50%. However, the recent CLEO 

n 

measurements |3| at Q^ = 13.48 GeV^, apparently with a very small error-bar of ±15%, can 
provide first possible opportunity to critically test theoretical predictions, although they 
do not shed light on the variation with Q^, which is a distinguishing feature of our result. 
Clearly, not only the moderate-energy time-like data seems reasonably reconciled, at higher 
energies both the CLEO result and the recent phenomenological prediction from J/tp decays: 



Mj/^ 



G'i^(M2= 9.59 GeV^ 



0.81 ± 0.18 GeV^ 



lie within reasonable range of our 



prediction for the total time-like form factor. This is surprisingly consistent with the pion 
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FIG. 7: Variation of the ratio of the pion and kaon time-hke form factors with Q^ in different 



approaches. The pion and kaon experimental data are taken from 



14. 



form factor results presented in 23|], also obtained within the light-cone /ct^- factorization 
approach, that agreed well with most of the available moderate-energy data (with statistics 



far better than the kaon data), including the CLEO data and a simi 



prediction 



46| based on J /ip decay analysis. Note that in the analysis 23| , the central value 



ar phenomenological 



of the twist-3 chiral parameter /Xtt was also taken to be 1.5 GeV. 

To this end, we consider the pion to kaon form factor ratios. In Fig. [TJ using the central 
result for the pion form factors from [23| (with /i^ = 1.5 GeV) we plot its variation with 



Q^ . The theoretical errors of the present work and 23|] are added in quadrature to obtain 
the error band as shown in the figure. The large error should not come as a surprise as the 
errors of both the pion and kaon factors are large. We now compare this result with other 
theoretical predictions and available experimental data. Note that the standard asymptotic 
pQCD result of Farrar and Jackson [7] yields a Q^ independent ratio. 



{F,G}r(Q^ 



^ = 0.67. 
Ik 



(40) 



Clearly, the central value of our time-like ratio deviates appreciably from the asymptotic 
value at low- and moderate-Q^, but however, it gradually approaches the asymptotic value 
at large Q^, and so does the space-like ratio. While our prediction fails to agree with the 
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very low-energy time-like data points [2| , showing the limitations of pQCD at such low-Q^ 
values, the higher Q^ ~ 4 GeV^ data points can somewhat be accommodated within our 
error-bars. At the same time, our time-like ratio at Q^ = 13.48 GeV^, i.e., 

G^(13.48GeV2) 



1.06 ±0.46, 



(41) 



Gi^(13.48GeV2) 

is surprisingly close to the CLEO value: 1.19 ± 0.17 at Q^ = 13.48 GeV^ [S], and the 
result obtained by taking the ratio of the phenomenologically estimated time-like pion form 
factor 46| and the time- like kaon form factor [39j], with the respective errors again added in 
quadrature: 



G.(M2 



J/i' 



9.48 GeV" 



Gk{M^ 



j/i, 



9.48 GeV' 



1.16 ±0.55 



(42) 



It is also noteworthy mentioning that the recent analysis [47[ based on a Light-Front Co- 
variant Model (LFCM) up to — g^ = 10 GeV^, yielded the ratio of the form factors quite 
similar to what we obtain in the space-like region. Finally, in Fig. [TJ we compare our re- 
sult, evidently working better towards large-Q^ values, with the soft QCD results obtained 
from QCDSR which are instead known to yield reliable predictions at low- and moderate-Q^ 
values. For example, the plot corresponding to the LD result [81] not only agrees well with 
the very low-Q^ space-like data (not resolved in the figure), but also with the low-energy 
time-like data when naively used in the time-like region. While, the analytically continued 



10| at low-energies yields a plot very different from that 



time-like LD result (see, Eq. 
of 8|, but toward larger-Q^ values both yield very similar predictions. Nevertheless, the 
QCDSR results significantly differ from the CLEO result and the one obtained from the 
phenomenological J/ip decay analysis. It is to be noted that in spite of the additional inclu- 
sion of the hard contributions, our space-like ratio of the total form factors does not differ 
significantly from that of |8|], except at the very low-Q^ ~ 0.2 GeV^ below which our result 
rapidly blows up. 

To sum up, in this paper we tried to systematically study the higher twist effects, namely, 
the twist-3 and twist-4 corrections to the standard twist-2 pQCD charged kaon form factors 
by adopting minimal model dependence arising from the inclusion of (a) the transverse 
degrees of freedom in the kaon wavefunctions/DAs, and (b) the non-factorizable soft QCD 
corrections via local duality. The work presented here extends and completes the analyses 
of the previous work 17|, |23|. Assuming the validity of the /c^- factorization ansatz through 
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the explicit TMD of the scattering kernel, we showed a non-trivial twist-3 contribution in 
the 2-particle sector which along with the large soft QCD corrections turn out to be the real 
hallmark of the "modified pQCD + soft QCD" approach to determine the space- and time- 
like kaon form factors. Other correction such as the 2-particle twist-4 were explicitly shows 
to have minor contributions only. To this end, the available moderate-energy experimental 
kaon data seems to be reasonably reconciled with the range of our predictions. It is also 
reassuring that the same approach works equally well independently for the electromagnetic 
pion form factors, which adds confidence to the arguments used in obtaining our results. 
It may, therefore, be speculated why the factorized result works so well for both the pion 
and kaon form factors in obtaining estimates, at least in the correct "ball-park" , in spite of 
factors like the resonances, hadronization and other final state interaction, naively neglected 
in this approach, that may render the factorized pQCD result questionable at the presently 
probed phenomenological region. However, to draw definite conclusion it is invaluable to 
have more high precision intermediate energy data, rather than to base our conclusions on 
such poor quality data. 

Acknowledgments: The authors would like to thank J.-W Chen and H.-N. Li and N.G. 
Stefanis for various discussions. 
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IV. APPENDIX 



2-particle Collinear Distribution Amplitudes (DAs) 



The 2-particle twist-3 collinear DAs for the charged kaon (say, K ) are defined at the 
scale of /io ~ 1 GeV, in terms of the following non-local matrix elements 25l427l| : 



\u[z)i-i^s[-z^ 



Jo 



{0\u{z)a^pj5s{-z)\K-{P)) 



/iA' < 1 



2' 



(1) 



rriu + m, 
rriK 
»i 



x{P^zp-Ppz^) / rfxe^«(^^V3;i^(a;,/io) 





with /ix = m\/{mu + ms) and ^ = 2x — 1. Note that the gauge-link factors (Wilson-line) in 
the matrix elements are to be implicitly understood. The normalization conditions for the 
above twist-3 DAs are given by 

fK 



^0 



2v/2iVc 



(2) 



which have the following asymptotic forms: 



,p(as) 
, (T (as) 



X) 



X] 



L 



K 



5fK 



V^K 



x{l — x) . 



The explicit formulas for the non-asymptotic 2-particle twist-3 collinear DAs, 



a series expansion over conformal spins at next-to- leading order, are given by 26 1 



(3) 



expressed as 



^I;i^(a;,/^^) 



+ ( -3r^3i.(/i^)c^3i.(/i^) - ^P^(-"') - ^Pi(/^')4(/i') ) Cy\0 



with 



a (as) / 



03;/.(^,/)=03;r(^) 1 + 



5773k(/U^ 

3 



:'73K(/^^)w3i^(/W^ 



20 



pl{f^')-^pUf^')a^{f^')]cnO 



VSK 



hKj_ 
Ik Pk 



PK 



Pk 



(4) 
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the non-perturbative parameters f^K and u^k being defined through the following matrix 
elements of local twist-3 operators: 

(0|«(0) a^xl^itD^, g^G^x] s(O) - ^ td^ m(0) a^xl5 9sGax siO)\K-{P)) 

3 

= T- ihKPaPuPi, ^3K + 0(higher twist) , (5) 

where gs is the strong coupling and Ga/3 is the gluon field tensor. The LO scale dependence 
of various twist-3 parameters are given by 



..(0) 



P;,(/.^) = L-3^//3op^(^2). ^U 



(0) 



VSKif^^) 


= r^>^V3M); 


ifl-\c,^\c,^ 


(^SKifJ-^) 


= L-tl^^uj,^[^l), 


7i2^4c. + lc. 


afif^') 


= L-^"/*af(p^); 


7r^^C., 



af(p^) = L-^"/A,4(^2^. 7f = |c^, (6) 

where L = as{fi'^)/as{^l), Cp = {N^ — l)/2Nc and Ca = A^c- However, the strange quark 
being massive, there is operator mixing of the ones in Eq. ([H]) with those of twist-2 operators, 
so that the resulting LO RG equations give the following scale dependences: 

+ ^ (^55/36,0 _^lWo)[;^^^^/C](^2)^ 

+ ^ {L''/^'''^^ - L^e/^/^o) [fj, m, af ](/.^) . (7) 

The 2-particle twist-4 coUinear DAs modify the twist-2 axial matrix element and are 
given by 

/■I r 1 

Jo I 4 

+ '-fKm%^z^ J dx e^«(^^) B4;x(x, fil) , (8) 
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where B4;/^ = g4.K — <p2;K, with the normahzation conditions expressed as 

^4*^' = I' dx {A, g},,K{x, /.2) = ^-^ , (9) 

and the asymptotic forms namely, 



A (as)/ \ ^^JK 2/1 \2 

A^.R- X = X (1 — x) . 



.SW^5^. (10) 

Next we display the explicit forms of the non-asymptotic twist-4 collinear DAs at NLO in 
conformal spin 26j: 

2\ 3/x ^., ^^[16 , 24 ^2N , on„ /,.2N , 20^ , 2N 



i4,Kix,fi ) = ^ ^_ x(l-x) j — + — a^ (/i ) + 20r73;^(^ ) + Y^4i^(^ ) 



135 
X {2x^(10- 15x + 6x^)lnx + 2x^(10- 15x + 6x^) Inx + xx(2 + 13xx)} 

,2N _ /^ /i , A , ^^.^^,2\ , fin^ ^,2n , 20 . oA ^1/2, 



where x = 1 — x, and in the notation of [25|, S^ = mj^ri4K and e = 21/8 ijJ4k- Note that the 
additional factor A/a ~ 3.5 in the denominator of A4./^ is a contrast to the expression given 
in 26|, which is introduced to normalize the DA. The non-perturbative parameters rj^x and 
U4K are defined through the following matrix elements of local twist-4 operators: 

~ 1 

{0\u{0)'JaigsG^uS{0)\K'{P)) = --fKm\ri4K iPf,gua - Pugf,a) , 



4 
9 



(0|n(0) [iD^, ig.G^xhx s{0) - - id^ n(0) igsG.x^x s(0)|ir-(P)) 



= fKm%V4K(^'iK ( PfiPiy - jm\ g^y j + C(twist 5) , (12) 

where, G^y = \e^upfjGP" is the dual gluon field tensor. Taking into account the mixing with 
the operators of lower twists, the LO RG evolution of the twist-4 parameters are 

r/4x(/i^) = V^-^''^V4k{i^ + i (l - U^tl^^ ■ 7S = I C 



3-^ 



v(0) /«„ r ■, , O. CO) 5 



[riAK U.k\ (/i^) = ^^^-/* \^.K ^.k\ (/ig) ; 7lt = -aCA. (13) 
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Finally, we present the various Gegenbauer polynomials used in the above formulas: 

cT{o = \(j>e-i), 

Cl'\0 = \aU'-^). (14) 
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